Cooling of superfluid neutron stars with muons 
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We extend our modeling of cooling of superfluid neutron stars by includ- 
ing the production of muons in the core, in addition to neutrons, protons, 
and electrons. The results are confronted with observations of middle-aged 
isolated NSs. Muons have little effect on the hydrostatic structure of NSs, 
on the slow cooling of low-mass NSs (RX J0822-43 and PSR 1055-52 in our 
model) and on the rapid cooling of massive NSs. They affect, however, the 
moderately fast cooling of medium-mass NSs (IE 1207-52, RX J0002+62, 
PSR 0656+14, Vela, and Geminga) and shift appreciably the mass range 
of these NSs to lower masses, which is important for correct interpretation 
of the observations. Moreover, the effects of muons can accurately be re- 
produced by a simple renormalization of NS models with no muons in the 
NS cores. 
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1. Introduction 

Rapid progress in detecting thermal emission from isolated neutron stars 
(NSs) with a new generation of orbital and ground-based observatories (e.g., 
Refs. 0, |2|) stimulates active theoretical studies of cooling isolated NSs. It 
is well known that cooling history of NSs depends on physical properties 
of matter of supranuclear density in NS cores. These properties (compo- 
sition of matter, equation of state - EOS, critical temperatures of various 
superfluids, etc.) are largely unknown: they cannot be reproduced in labo- 
ratory or unambiguously calculated (because of the lack of exact manybody 
theory for describing the systems of particles interacting via strong forces) . 
However, they may be constrained by comparing cooling calculations with 
observations of isolated NSs (e.g., Ref. ||). 

We will focus on the theoretical interpretation of observations proposed 
recently in Refs. @, |, § (hereafter KYG), and g, §. The authors re- 
stricted themselves to model EOSs in the NS cores in which the presence 
of muons was neglected. In the present paper we show that the effect of 
muons on the cooling may be important. 

For the observational basis, we take the same isolated middle-aged NSs 
as KYG, but exclude RX J1856-3754. The age of this intriguing radio-quiet 
NS has recently been revised ||; the present value i = 5x 10 5 yrs is lower 
than the former one, 9 x 10 5 yrs, and the source has become less restrictive 
for cooling theories (cf. KYG and Ref. 0). In addition, there are indications 
that the emission from the entire NS surface is obscured by the emission 
from a hot spot on the surface of the NS; if so the surface temperature is 
poorly determined from the present observations (e.g., 0] and references 
therein) . 

The effective surface temperatures, T s °°, redshifted for a distant observer, 
and ages t of seven isolated NSs are taken from Table 3 of KYG and dis- 
played in Figs. 2 and 3 below. The three youngest objects, RX J0822-43 
|0|, IE 1207-52 and RX J0002+62 g§, are radio-quiet NSs in super- 



nova remnants. 

The other objects, Vela (PSR 0833-45) H, PSR 0656+14 pi, Geminga 



(PSR 0633+1748) @, and PSR 1055-52 JT^ are observed as radio pulsars. 
The adopted values of T s °° are inferred from the observed spectra using 
various models of stellar emission described in KYG. Recently, the values 



of T s °° for some of the sources have been revisited in Refs. pi 0]> |18|- p4|. 
Since the new data are basically in line with those used in KYG, we do not 
introduce the corresponding changes; they are insignificant for our analysis. 

As shown in KYG, the observations can be explained using the models of 
NSs with the cores composed only of neutrons, protons, and electrons, and 
assuming the presence of nucleon superfluidity with the density dependent 
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critical temperatures T c (p). Following KYG we consider superfluidities of 
three types produced by: (1) singlet-state pairing of protons in the NS 
core (T c = T cp ); (2) singlet-state pairing of free neutrons in the inner crust 
(T c = T cns ); and (3) triplet-state pairing of neutrons in the core (T c = T cnt ). 
Owing to a large scatter of microscopic theoretical models of T c (p) (e.g. 
Ref. p5|), we treat T c (p) as free parameters. 

KYG considered cooling of NSs with rather strong pairing of protons and 
weak pairing of neutrons in the core, and with a strong pairing of neutrons 
in the crust. They found that cooling middle-aged NSs can be divided into 
three types. 

• Type I NSs are low-mass NSs which show slow cooling with (modified 

or direct) Urea processes of neutrino emission strongly suppressed by 
proton superfluidity. The cooling curves, T s °°(i), are insensitive to NS 
mass, EOS in the core, and proton superfluidity (i.e., to T cp ) as long 
as the latter is sufficiently strong. KYG interpreted RX J0822-43, 
and PSR 1055-52 as low-mass NSs. 

• Type II NSs are medium-mass NSs which show moderately fast cooling 

regulated by direct Urea process partly reduced by proton superflu- 
idity in the NS central kernels. The cooling curves are sensitive to 
NS mass, EOS, and especially the T cp (p) profiles in the NS kernel. If 
the EOS and T cp (p) are fixed, the effective surface temperature de- 
creases smoothly with increasing M, and one can measure the mass 
('weigh" medium-mass NSs) using the observed limits on T s °°(i). KYG 
treated IE 1207-52, RX J0002+62, Vela, PSR 0656+14, and Geminga 
as medium-mass NSs. 

• Type III NSs are massive NSs which show rapid cooling via direct Urea 

process in the NS kernels, almost unaffected by proton superfluidity. 
The surface temperatures of these NSs are low (a few times 10 5 K for 
t ~ 10 3 yrs), being not too sensitive to the NS structure. No NS of 
such type has been observed so far. 



2. Cooling models 

We use the same cooling code as in KYG and modify the physics input 
in the NS core to include the effects of muons. 

First, we have included muons in the EOS. We use a stiff EOS proposed 
in Ref. 26], the model I of the symmetry energy of nucleon matter with 
the compression modulus of saturated nuclear matter K = 240 MeV. The 
same model of nucleon-nucleon interaction was adopted by KYG (EOS A in 
their notations) who, however, artificially suppressed the presence of muons. 
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Table 1. NS models employing EOSs A and Au (without and with muons) 



Model 


Main parameters 


EOS A 


EOS Ap 


Maximum 




1.977 


1.950 


mass 


Pcmax/10 14 g cm" 3 


25.75 


26.55 


model 


R, km 


10.754 


10.602 


Direct Urea 


M D /M 


1.358 


1.249 


threshold 


Pd/10 14 g cm -3 


7.851 


7.423 


model 


R, km 


12.98 


12.952 



Now we include the muons and obtain EOS A//. We will compare the results 
obtained with EOSs A and Ap,. 

For EOS Ap,, the muons appear at p > p^ = 2.55 x 10 14 g cm -3 (when the 
electron chemical potential exceeds the muon rest energy). Their fraction 
is lower than 10% everywhere in the NS core. Their appearance slightly 
softens the EOS, slightly increases the fraction of protons and decreases the 
fraction of electrons. These changes weakly affect the NS internal structure, 
reducing the maximum NS mass by 1.4%. The masses M, central densities 
p c , and radii R of two NS configurations for EOSs A and Ap are given in 
Table 1. The first configuration corresponds to a maximum-mass NS. The 
second corresponds to the onset of powerful direct Urea neutrino emission 
[§7| in the NS kernel (p c = p D , M = M D ) . 

Another modification of the cooling code consists of incorporating the 
heat capacity of muons. It is straightforward (e.g., Ref. |2^]) and has no 
noticeable effect on the NS cooling. 

Finally, we have modified the emissivities of neutrino reactions in the NS 
cores (as described in Ref. (^9|). The main effect of muons is to lower the 
threshold density pd of onset of direct Urea process (Table 1). The lowering 
is associated with the changes of the fractions of neutrons, protons, and 
electrons in muonic matter which relaxes the onset condition. At p < p^, 
modified Urea process is basically the leading one. At p > p^, two new 
branches of this process appear, the neutron and proton branches, where 
the muons participate instead of the electrons. However, these new branches 
have little effect on the cooling because their emissivity is comparable to 
the emissivity of the ordinary modified Urea process. The ordinary direct 
Urea process operates at p > p-Q. At still higher density p > pd^ (pDfi = 
9.257 x 10 14 g enr 3 , for EOS Ap) another direct Urea process is open, where 
the muons participate instead of the electrons. We have included it into the 
code but, again, it has little effect on the cooling because its emissivity is 
comparable to the emissivity of the ordinary direct Urea process. 

Now we describe our models of superfluidity in NSs. As in KYG, we 
assume that the triplet-state pairing of neutrons in the NS cores is weak 
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Fig. 1. Density dependence of the critical temperature of proton superfluidity 
for models lp and 2p, and for EOS Ap (with muons, thick lines) and EOS A 
(without muons, thin lines). Vertical dotted lines indicate the direct Urea (Durca) 
threshold. 

(maximum T cn t < 10 s K) and can be ignored. Otherwise, a strong neutrino 
emission due to Cooper pairing of neutrons would accelerate the cooling 
and hamper the interpretation of observations of older NSs. The singlet- 
state pairing of neutrons in the crust controls the neutrino luminosity of 
low-mass NSs and has no direct relation to the presence of muons in the NS 
core. While exploring the effects of muons we ignore this pairing (see also 
Sect. 3). 

Thus, we focus on the proton superfluidity in the core. We parameterize 
the dependence of T cp on the proton Fermi wavenumber k = (37r 2 n p ) 1//3 
(ra p being the proton number density) by the same expression as Eq. (1) 
in KYG and consider two models of superfluidity, lp and 2p, employed by 
KYG. The corresponding T cp (p) profiles are plotted in Fig. |l|. The maximum 
of T cp (p) is about the same (~ 7 x 10 9 K) for all models, but superfluidity 
2p extends to a higher densities than lp. Since we use the same model of 
nucleon interaction for EOSs A and Ap, we have the same dependence of 
T cp on k (but different dependences of T cp on p because the proton fractions 
are different for EOSs A and Ap). The curves in Fig. 1 are typical of the 
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Fig. 2. Observational limits on surface temperatures of seven NSs compared with 
cooling curves for NSs of several masses, adopting lp proton superfluidity and either 
EOS A/i (with muons, solid lines) or ESO A (without muons, dashed lines). If two 
masses are given, an upper value corresponds to EOS A and a lower (boldface) 
value to EOS Ap. 

microscopic calculations (e.g., Ref. p5| ). As discussed above, the muons 
shift both the direct Urea threshold and the T cp (p) profile to lower values 
of p. The latter shift, due to the increase in the proton fraction, is bigger. 
For instance, T cp (p) for model lp of proton superfluidity vanishes below the 
direct Urea threshold in the muonic matter, but persists above the direct 
Urea threshold in the matter without muons. 

3. Results and discussion 

The results of the cooling calculations are illustrated in Figs. || and ||, 
which also show the observational data. 

According to our models, the cooling curves of low- mass NSs with and 
without muons are almost indistinguishable. The neutron superfluidity in 
the crust can affect the cooling of these stars in exactly the same way as 
discussed in Ref. || and KYG for the NSs without muons. The calculations 
also show that muons have almost no effect on rapid cooling of massive NSs. 
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Fig. 3. Same as Fig. 2, but for proton superfluidity model 2p. 

However the effect of muons on moderately fast cooling of medium-mass 
NSs is substantial. Figure 2 corresponds to model lp of proton superfluidity. 
Without muons, the T cp {p) profile extends above the direct Urea threshold 
(Fig. 1). Accordingly, we have a representative class of medium-mass NSs 
(dashed curves). By varying the mass M from 1.36 M to 1.465 M we 
can explain the observations of the five sources (IE 1207-52, RX J0002+62, 
Vela, PSR 0656+14, and Geminga) by the cooling of medium-mass NSs, and 
thus 'weigh" these NSs (Ref. [fl, KYG). In the presence of muons, the proton 
critical temperature T cp {p) drops to zero below the direct Urea threshold, 
Pd (Fig- l)i he., the proton superfluidity disappears in the centers of NSs 
with masses somewhat lower than the threshold mass Md- The transition 
from slow to fast cooling occurs in a very narrow mass range, from 1.249 M© 
to 1.264 Mq, so we have no representative class of medium-mass NSs. The 
proposed interpretation of the five sources as medium-mass NSs containing 
muons becomes unlikely (just as in the absence of superfluidity, see KYG). 

Figure 3 is the same as Fig. 2 but for model 2p of proton superfluid- 
ity. This superfluidity extends above p-Q with and without muons (Fig. 1). 
Therefore, we have representative classes of medium-mass NSs and we can 
weigh our five medium-mass sources in both cases. Without muons, as in 
KYG, we obtain the masses from 1.49 Mq to 1.60 M . With muons, we 
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obtain significantly lower (and narrower) mass range for the same sources, 
from 1.27 M to 1.34 M . 

The inferred mass range changes in the presence of muons because they 
lower the threshold density of direct Urea process and shift the decreasing 
slope of the T cp (p) profile to an even lower density (Fig. 1). Accordingly, 
direct Urea process becomes stronger at lower p and this accelerates the 
cooling of lower-mass NSs. Therefore, the mass range of medium-mass NSs 
shifts to lower M. 

4. Conclusions 

We have simulated the cooling of superfluid NSs with and without 
muons. The inclusion of muons does not violate the main results of Refs. 
(§]-|| on the existence of three types of cooling superfluid NSs. As in KYG, 
we can interpret RX J0822-43, and PSR 1055-52, as low-mass NSs, and IE 
1207-52, RX J0002+62, Vela, PSR 0656+14, and Geminga, as medium- 
mass NSs. This interpretation requires strong singlet-state proton pairing 
and weak triplet-state neutron pairing in the NS cores. 

Muons do not change the hydrostatic structure of NSs. The muon heat 
capacity and direct/modified Urea processes involving muons have almost 
no effect on NS cooling. The cooling curves of low-mass and high-mass NSs 
are insensitive to the presence of muons. 

However, muons do affect the cooling of medium-mass NSs and the 
associated 'weighing" of observed NSs proposed in Ref. Q and KYG. For 
the same nucleon-nucleon interaction model and the superfluid properties 
of NS cores, we obtain noticeably lower masses of NSs with muons than 
without muons (Fig. 3). Therefore, one should take muons into account for 
a correct interpretation of observations. 

It is remarkable that muons affect the cooling by their mere existence, 
as passive particles. We could switch off all neutrino reactions and heat 
capacity associated with muons (i.e., retain the 'old" physics input in the 
NS core) but employ the fractions of neutrons, protons and electrons in the 
EOS modified by the presence of muons. The cooling curves obtained in 
this numerical experiment would almost coincide with the exact ones. This 
means that the effects of muons can be accurately incorporated by renormal- 
izing the NS models containing no muons at all. In other words, any model 
of a cooling NS with muons is equivalent to a certain model of a cooling NS 
without muons. Clearly, the muons will have similar effects in the presence 
of hyperons in supranuclear matter. This circumstance may be helpful in 
constructing the codes to simulate cooling of NSs containing hyperons. 

The existence of three types of cooling NSs is effectively regulated by two 
factors. The first is the onset of the powerful direct Urea process of neutrino 
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emission at densities p > pr>', the second is the strong proton superfluidity 
extending at p > pr>. The latter suppresses the direct and modified Urea 
processes at or just above the critical density p-& but opens them at higher 
p. In this way the superfluidity smears out the sharp transition between 
the slow and rapid cooling regimes and creates a representative class of 
medium-mass NSs, favorable for the interpretation of observations. These 
conditions may be also realized in NS cores containing hyperons, pion or 
kaon condensates, or quark matter. We will analyze different models of NS 
cores in future publications. 

We are grateful to G.G. Pavlov for encouragement, to D. Page for criti- 
cism, and to P. Haensel and A.D. Kaminker for careful reading of the text. 
The work was partially supported by RFBR (grant No. 00-07-90183) and 
KBN (grant No. 5P03D.020.20). 
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